int niceGraphs(){
TFile *f = new TFile("sbruResults.root"); 
gStyle->SetOptStat(kFALSE);
TCanvas* c2 = new TCanvas("c2");
TH2F *h2 = (TH2F*)f->Get("run_4/monitorShell_g_en_theta_3;1");
h2->GetXaxis()->SetTitle("Energy(MeV)");
h2->GetXaxis()->SetRangeUser(0,2.1);
h2->GetZaxis()->SetRangeUser(0,10);
h2->GetYaxis()->SetTitle("Angle(deg)");

//h2->GetZaxis()->SetDrawOption("")
h2->SetTitle("");
h2->Draw("COLZ");

TLine *line = new TLine(0.0,85.0,2.1,85);
   line->SetLineWidth(2.2);
   line->Draw();

TLine *line = new TLine(0.0,95.0,2.1,95);
   line->SetLineWidth(2.2);
   line->Draw();

TLine *line = new TLine(0.0,130.0,2.1,130);
   line->SetLineWidth(2.2);
   line->Draw();

TLine *line = new TLine(0.0,140.0,2.1,140);
   line->SetLineWidth(2.2);
   line->Draw();

c2->Modified();
};
int niceGraphs2(){
	TFile *f= new TFile("5MeV.root");
gStyle->SetOptStat(kFALSE);

TCanvas* c2 = new TCanvas("c2");
c2->Divide(2,1);
c2->cd(1);
TH1F *h1 = (TH1F*)f->Get("ThetaSign2;1");
h1->GetXaxis()->SetTitle("Angle (rad)");
h1->GetYaxis()->SetTitleOffset(1.5);
h1->GetYaxis()->SetNdivisions(5,kTRUE);
h1->GetYaxis()->SetTitle("Number of scattered photons");
h1->GetXaxis()->SetRangeUser(-0.2,0.2);
//h1->GetYaxis()->Set
c2_1->SetLogy();
h1->GetXaxis()->SetRangeUser(-1,1);
h1->SetTitle("");
h1->Draw();
h1->GetXaxis()->SetRangeUser(-0.02,0.02);
c2->Update();


c2->cd(2);
TH1F *h2 = (TH1F*)f->Get(" E(eV);1");
h2->GetXaxis()->SetTitle("Photon energy (keV)");

h2->GetYaxis()->SetTitleOffset(1.5);
h2->GetYaxis()->SetNdivisions(5,kTRUE);

h2->GetYaxis()->SetTitle("");
h2->SetTitle("");
h2->Draw();
h2->GetXaxis()->SetRangeUser(0,5100);
c2->Update();

//h1->GetZaxis()->SetRangeUser(0,10);
//h1->GetYaxis()->SetTitle("Angle(deg)");
//h1->Draw();
};


